Generalizations of Wiener polarity index and 
terminal Wiener index 



Aleksandar Ilic * 
Faculty of Sciences and Mathematics, Visegradska 33, 18 000 Nis 
University of Nis, Serbia 
e-mail: aleksandari@gmail . com 

Milovan Ilic 

Faculty of Information Technology, Trg republike 3, 11 000 Beograd 
University of Belgrade, Serbia 
e-mail: ilic.milovanOgmail.com 

June 16, 2011 
Abstract 

In theoretical chemistry, distance-based molecular structure descriptors are used for 
modeling physical, pharmacologic, biological and other properties of chemical compounds. 
We introduce a generalized Wiener polarity index Wk (G) as the number of unordered pairs 
of vertices {u, v} of G such that the shortest distance d(u, v) between u and v is k. For 
k = 3, we get standard Wiener polarity index. Furthermore, we generalize the terminal 
Wiener index TWk(G) as the sum of distances between all pairs of vertices of degree k. 
For k = 1, we get standard terminal Wiener index. In this paper we describe a linear 
time algorithm for computing these indices for trees and partial cubes, and characterize 
extremal trees maximizing the generalized Wiener polarity index and generalized terminal 
Wiener index among all trees of given order n. 

Key words: Distance in graphs; Wiener polarity index; terminal Wiener index; Wiener 
index; Partial cube; Graph algorithm. 
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1 Introduction 

Let G = (V,E) be a connected simple graph with n = \V\ vertices and m = \E\ edges. For 
vertices u, v G V, the distance d(u, v) is defined as the length of the shortest path between u 
and v in G. The diameter diam{G) is the greatest distance between two vertices of G. Let 
dk(u) denotes the number of vertices on distance k from the vertex u. Let deg(v) denotes the 
degree of the vertex v. 

In theoretical chemistry molecular structure descriptors (also called topological indices) 
are used for modeling physico-chemical, pharmacologic, toxicologic, biological and other prop- 
erties of chemical compounds |12| . There exist several types of such indices, especially those 
based on vertex and edge distances. Arguably the best known of these indices is the Wiener 
index W, defined as the sum of distances between all pairs of vertices of the molecular 
graph [9] 

W(G)= Yl d ^ v )- 

u,v€V(G) 
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Besides of use in chemistry, it was independently studied due to its relevance in social science, 
architecture, and graph theory. With considerable success in chemical graph theory, various 
extensions and generalizations of the Wiener index are recently put forward p2 \TT\ 124"] . 

The Wiener polarity index of a graph G is defined as the number of unordered pairs of 
vertices {u, v} of G such that the shortest distance d(u, v) between u and v is 3, 

WP{G) = \{{u,v) | d(u,v) =3, u,ve V}\. 

Hosoya [13] found a physico-chemical interpretation of WP. Du, Li and Shi |8j described a 
linear time algorithm for computing the Wiener polarity index of trees, and characterized the 
trees maximizing the index among all trees of given order. Deng et al. |U [5j [6] and Liu et 
al. [21] characterized extremal n-vertex trees with given diameter, number of pendent vertices 
or maximum vertex degree. 

The terminal Wiener index of a graph G is defined by Gutman, Furtula and Petrovic in 
[TO] as the sum of distances between all pairs of pendent vertices of G 

TW(G)= Y, 

deg(u)=deg(v)=l 

Furthermore, the authors described a simple method for computing TW of trees and char- 
acterized the trees with minimum and maximum TW. Heydari and Gutman [13] provided 
a formula for calculating TW of thorn graphs, while Deng and Zhang [7] studied equisepa- 
rability on terminal Wiener index. Independently, Szekely et al. [23 introduced the same 
index (the sum of distances between the leaves of a tree) and studied the correlation between 
various distance-based topological indices. 

The paper is organized as follows. In Section 2 we introduce generalization of the Wiener 
polarity index and characterize the trees maximizing the generalized Wiener polarity index 
among all trees of given order, while in Section 3 we designed linear algorithm for calculating 
this index. In Section 4 we introduce generalization of the terminal Wiener index and charac- 
terize trees maximizing the generalized terminal Wiener index among all trees of given order 
and k < 3. In Section 5 we present formula for calculation of TWk(G) for partial cubes and 
in particular closed formula for TW3 of coronene series H^. We close the paper in Section 6 
by proposing new problems for research. 



2 Generalization of Wiener polarity index 

For k > 1, we define the generalized Wiener polarity index as the number of unordered pairs 
of vertices {u, v} of G such that the shortest distance d(u, v) between u and v is k, 

veV(G) 

Notice that W{G) = Y,t=? {G) W k {G). 

For k = 1, it can be easily seen that W\{G) = m, where m is the number of edges. For 
k = 2, we have 

W 2{G) = £ = Sgj^W _ m _ Ml(0) _ m , 

where M\{G) denotes the first Zagreb index of a graph |22j. 
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For k = 3 we have the Wiener polarity index, 

W ^ T ) = ^ de 9( v ) ~ l )(deg(u) - 1) = deg(u)deg(v) - deg 2 v + m 

uv£E uv£E vEV 

= M 2 (T)-M 1 (T) + m, 

where M 2 (T) denotes the second Zagreb index of a graph |15j . 

In the following assume that k > 3. If the diameter of T is less than k, then W k (T) = 0. 
Therefore, the minimum value of W&(T) is zero, and it is achieved for all trees with diam(T) < 
k (for example the star S n ). On the other hand, we will prove that the maximum value of 
Wfc(T) is achieved for a tree with diameter k and with all pendent vertices on distance k. 

The group of pendent vertices is defined as the set of all pendent vertices attached to the 
same unique neighbor. Let Ai and A 2 be two different groups of pendent vertices with the 
unique neighbors w\ and w 2 , such that the distance between arbitrary pendent vertices from 
these groups is not equal to k. Let 

Pi = d k {v) and p 2 = d k (u). 
veAi u&A 2 

Without loss of generality assume that p\ > p 2 . If we remove all pendent vertices from A 2 
and add them to the set A\, we get a new tree T' such that 

W k (T')-W k (T)= Pl -p 2 >0. 

By repetitive application of this transformation, we will get a new tree with possibly increased 
generalized Wiener polarity index. The diameter of T' is not greater than the diameter of T' 
and each transformation introduces at most one new pendent vertex. By choosing two most 
distant groups of pendent vertices, we will get the extremal tree with diameter equal to k 
and all pendent vertices are on distance k or 2. 

Assume that there are p groups of pendent vertices with sizes a\ , a 2 , ■ ■ ■ , a p and a\ + a 2 + 
. . . + a p = q. Since diam{T) = k, we have n — k + 1 > q > 2. The distance between any two 
pendent vertices not from the same group is equal to k, and therefore 

8=1 V 1=1 

The minimum value of X^f=i a f un der the condition Yli=i o-i = q is achieved if and only if all 
numbers aj are as close as possible, i. e. |aj — aj\ < 1 for all 1 < i < j < p. This can be 
easily proved by the transformation (aj, a,-) i— >• (aj + 1, cij — 1) with a,- > aj + 2, since 

(at + l) 2 + (aj - l) 2 - a 2 - a 2 = 2(oj - o 3 -) + 2 > 0. 

Notice that the tree is uniquely determined by the distances between pendent vertices 
|26j . A starlike tree is a tree with exactly one vertex of degree at least 3. If p = 2, we have 
W k (T) = a\a 2 and a\ + a 2 = n — k + 1 and finally 

n — k + 1 
2 

Let p > 2. For k odd, we can consider two groups of pendent vertices together with the 
unique path connecting them. The third group of pendent vertices must be on equal distance 



W k (T) 



n 



k + l 
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from both groups and that is impossible. Therefore, the extremal value for odd k is achieved 
for p = 2. 

For k even, similarly it can be concluded that there is a unique tree with p groups 
of pendent vertices (starlike tree with p paths with equal lengths together with groups of 
pendent vertices attached at the pendent vertices of these p paths). For p > 2, we have 

n=l+pQ-l) +Y^Oi = l+p^-lj +q, 

and since p < q, we have p < 2 • Therefore, using Cauchy-Schwartz inequality it follows 



W k (T) = ^(V-X> 2 ) 



< 7J 9 

2 V V 



l L-i- p ±+p 



i - 



p 



Let 



,(,)_» („-!-!*+,) (l-i 



for 2 < _p < 2 • ^ < 2 • . The first derivative equals 

(2 - 2n - 2p + fep) (2 - 2n + 2p - kp - 4p 2 + 2kp 2 ) 



fip) 



4p2 



Finally it holds that f(p) is increasing for 2 < p < p* and decreasing for p* < p < 2 ■ ^r-, 
where 



p 



k-2 + ^J(k - 2)(16n + k - 18) _ 1 1 ll&n + k- 18 

~ 4 + 4 V Jfc-2 



4(fc - 2) 

is the second largest root of f'(p) = 0. Therefore, the maximum of the Wiener polarity index 
for k even is achieved for p = 2 or p that is near to p* (see Figure 1). 




Figure 1: The extremal tree with maximal generalized Wiener polarity index on n = 22 and 
k = 6. 
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3 Linear algorithm for Wk of trees 



Let T be an arbitrary tree rooted at the vertex 1. We will process the vertices according 
to the distance from the root vertex or in the recursive depth first search method [2J. For 
each vertex v, we will keep the vector a[v] of the length k + 1 that stores the number of 
vertices in the subtree under v on distances 0,1,2, ... ,k. It follows that a[i>][0] = 1 and 
a[f][l] = deg(v) — 1 for all vertices different then the root. The matrix a with dimensions 
n x k is computed recursively in the first procedure. 

In the second procedure we calculate the generalized Wiener polarity index. In DFS tree, 
for each path of the length k there is a unique vertex v on the smallest distance from the 
root. Therefore, we need to traverse all vertices v and count the vertices that are on distance 
k in the subtree under v, such that the unique path connecting these vertices contains v. In 
the first case, we just add 2a [u] [A;]. In the second case, we need to consider all neighbors u of 
v different than parent[v] and for each j = 0, 1, . . . , k — 2 count the number of vertex pairs 
(x, y) such that: 

• x is in the subtree under u and y is not; 

• x is on the distance j from u; 

• y is on the distance k — j — 1 from v and under v. 

Finally we counted every vertex pair on distance k twice, as showed in the second procedure. 



Procedure DFS (vertex v) 
Input: The adjacency list of the tree T with the root vertex root . 
Output: The array parent and matrix a. 

a[v][0] = l; 

foreach neighbor u of v do 

if (parent[u] = 0) and (u ^ root) then 
parent[u] = v; 
DFS(u); 

for i = to k — 1 do 

I «[«][* + 1] =o[«][i + i] + <*[«][*]; 

end 
end 
end 



4 Generalization of terminal Wiener index 

For k > 1, we define the generalized terminal Wiener index as the sum of the distances 
between all pairs of its vertices of degrees k, 

TW k {G)= d M- 

u,v£V(G) 
deg(u)=deg(v)=k 

For k regular graphs, WTk{G) = W(G) and WTi(G) = for i ^ k. For k = 1, we have 
terminal Wiener index. 
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Procedure GWP 

Input: The adjacency list of the tree T with the root vertex root and the array a . 
Output: The generalized Wiener polarity index Wk. 



for v = 1 to n do 

Wk = Wk + 2-a[v][k]; 
foreach neighbor u of v do 
if parent [v] ^ u then 
for i = to k — 2 do 
[ Wk = Wk + a[u][{\ ■ (a[v][k - 1 - i] - a[u][k - 2 - *]); 
end 
end 
end 
end 

return Wk/2; 



Theorem 4.1 (|9j) Let T be a tree on n vertices. Then, 



with equality if and only ifT = S n or T = P n . 

Theorem 4.2 ( |16| ) Let w be a vertex of a nontrivial connected graph G. For nonnegative 
integers p and q, let G(p, q) denote the graph obtained from G by attaching to vertex w pendent 
paths P = wv\V2 ■ ■ ■ v p and Q = wu\U2 ■ ■ ■ u q of lengths p and q, respectively. If p > q > 1, 
then 



For k = 1, Gutman, Furtula and Petrovic in [TO] characterized the extremal trees within 
the class of all trees with n vertices that maximize terminal Wiener index. For k = 2, the 
unique tree that maximizes TW2 is the path P n . Interestingly, the path P n is the unique tree 
that minimizes TW\ . This follows from Theorems 14.11 and 14.21 and the simple fact that in 
every tree there are at least 2 pendent vertices. 

Let T* be the extremal tree that maximizes the generalized terminal Wiener index for 
k > 3. It can be easily proved that there are no pendent paths P = v if 2 . . . v p of length 
greater than 2 attached at some vertex of T*. Otherwise, remove pendent edges one by one 
and subdivide any edge e such that both components of T* — e contain vertices of degree k - 
this way we increase TWk- If such edge does not exist, remove the edge v v -\v v and add new 
edge v p -2Vp - this way the generalized terminal Wiener index remains the same, or increases 
if and only if k = 3. 

Lemma 4.3 Among trees on n vertices the maximal number of vertices of degree k > 2 is 



Proof. Consider the induced graph H composed of the vertices of degree k. Let h be the 
number of vertices in H, and let / be the number of edges in H. Since H is acyclic and 



Wk = 0; 




W(G(p,q))<W(G(p + l,?-!)). 



m(n, k) 



n-2 
k-1 
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possible disconnected, we have / < h — 1. Furthermore, for the number of edges in tree T 
holds 

h-k-(h-l)<h-k-f<n-l. 
It follows that h{k — 1) < n — 2, which completes the proof. ■ 

v vv , , vv 

Figure 2: The extremal tree with maximal generalized terminal Wiener index on n = 20 and 
k = 4. 

The C nt k,p be the caterpillar obtained from a path of length s + 2 = n — p ■ (k — 2), by 
attaching k — 2 pendent vertices to exactly p vertices of a path P s +2 = ^o^i • • • v s v s+ i, starting 
from the vertices v\ and v s on the both ends towards center symmetrically (see Figure 2). 
The generalized terminal Wiener index can be easily calculated, 

TW k (C n: k,p) = ^2 d(vi,Vj) 

l<i<[p/2] 
s-[p/2\+l<j<s 

= (p - - 1) + (p - 3)(a - 3) + (p - 5)(s - 5) + . . . 
J ^p(3ps — p 2 — 2), if p is even 

I T2~(P + -OOP — l)(^ s — p)> if p is odd. 

/ (3 n P + 5p 2 — 3A;p 2 — 2 — 6p) , if p is even . . 

1 12 (P — l)(3ra + 5p — 3A;p — 6), if p is odd. 




Figure 3: Transformation that increases generalized terminal Wiener index TW%. 

We call a caterpillar C 3-bounded if all vertices of C have degree less than or equal to 3. 
Theorem 4.4 Let T be a tree on n > 4 vertices. Then 

TW 3 (T) < TW 3 (C nMn/2hl ), 
with equality if and only if T = C n 3 i n /2j-i- 

Proof. Let T* be a rooted tree with maximal value of generalized terminal Wiener index for 
k = 3. If T* is not a caterpillar, consider a branching vertex to such that in the subtree under 
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w there are only 3-bounded caterpillars attached at w. Let C\ and C2 be two caterpillars 
attached at w, such that C\ has p vertices v±,V2, ■ ■ ■ ,v p of degree 3 and C2 has q vertices 
ui, U2, ■ ■ ■ , u q of degree 3 (see Figure 3). Without loss of generality, we can assume that the 
number of vertices of degree 3 in G = T* \ {61,62} is greater than or equal to p and q, 
namely r > p > q. 

Let ai,a2, ■ ■ ■ a p be the distances from the vertex w to the vertices v±, V2, ■ ■ ■ , v p and 
61, 62, • • • , b q be the distance from the vertex w to the vertices u\, U2, ■ ■ ■ ,u q . Let D(w) be 
the sum of distances from the vertex w to the vertices of degree 3 in the subgraph G. 

The generalized terminal Wiener index of the tree T* equals 



TW 3 (T*) = TW z {G) + {p + q)D{w)+r\Y,a l + Y J h 

i=i 3=1 



+ ^2( a J ~ ai ) + ^2( b J -bi) + q-^2ai+p-^2b 

i<j i<j i=l j=l 



After reattaching the caterpillar C2 to the end of caterpillar C\, the degree of vertex w 
remains the same. Furthermore, the vertex v p has degree 3 and the generalized terminal 
Wiener index of transformed tree T' equals 



TW 3 (T') = TW 3 {G) + {p + q)D(w) + r-^2a i + riq-a p + J2 b 

t=l I 3=1 



q 

^3 



+ ^( a j ~ a ») + ^2( b J -h) + q- I V ■ a P ~ ^2 ^ j + p ■ ^ b 

i<j i<j \ i=l / j=l 



p \ q 



By subtraction, we get 

p 

TW 3 (T') - TW 3 (T*) = rq -a p +pq-a p -2q-Y, <H, 

i=l 

Since r > p and a p > ai for i = 1, 2, . . . ,p, we have 

TW 3 (T')-TW 3 (T*) = gfo p (r+p)-2j3o i J 



> 2q ■ Y^{a p - a,) > 0, 



i=l 

with equality if and only if p = q = 1. 

It follows that using this transformation we can not decrease the generalized Wiener index 
TW 3 . If deg(w) > 3, we can remove the pendent edge and subdivide some edge e such that 
both components of T* — e contain vertices of degree 3, and further increase the generalized 
Wiener index. Finally, we conclude that the extremal tree is a caterpillar of the form C n ^ )P - 

For k = 3, from ([!]) we have 



and 



J2P (3rap — 4p 2 — 2 — 6p) , if p is even 

j2-(p + l)(p — l)(3n — 4p — 6), if p is odd. 



fip) ~ Sip - 2) 



(p — l)(n — 2p) — 1, if p is even 

{p — l)(n — 2p), if p is odd, 
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which is greater than zero since n > 2p + 2. By direct verification, we get /(|_§J — 1) > 
/(LfJ -2) forn>4. 
Therefore, it holds 

W 3 (T)<W3(C n ,3, W2j -l), 
with equality if and only if T = C n>s i n /2j-i- ^ 



5 Calculating TWk of partial cubes 

Let G be a connected graph. Then e = xy and f = uv are in the Djokovic- Winkler relation 
[25] if 

<i(x, u) + d(y, v) ^ d(x, v) + d(y, u). 

The relation is always reflexive and symmetric, and is transitive on partial cubes. 
Therefore, partitions the edge set of a partial cube G into equivalence classes F\, F%, . . . , F s , 
called G-classes. For any 1 < % < s, the graph G — Fi consists of two connected components. 
The vertex sets of these components will be denoted with Wio and Wn, because they can 
described as the vertices whose i-th coordinate is and 1, respectively. The sets W{ tX , 
1 < i < s, x £ {0) 1} are called halfspaces of G, while W^o and Wi t \ are complementary 
halfspaces [18j and it holds \W {ifl) \^ + |W (iil) \ {k) = 1. 

Let G be a partial cube with halfspaces Wi iX , 1 < i < s, x £ {0, 1}. For any 1 < i < s 
and any x € {0, 1}, I ^ / (i s x) I denotes the number of vertices with degree k. 

Theorem 5.1 Let G be a partial cube with halfspaces W^ x , 1 < % < s, % € {0, 1}. Then 



T^(G) = ^|^ i0) |W-|^ 1) |( fe ). 



i=l 



Proof. G is a partial cube, hence vertices of G can be considered as a binary s-tuple u = 
u\U2 ■ ■ ■ u s . Moreover, since G is isometric in Q s , the distance between two vertices is the 
number of positions in which they differ. Set 5(x, y) = if x = y, and 5(x, y) = 1 for x ^ y. 
Then 



TW k (G) 



E 



d(u, v) 



u,v€V(G) 
deg(u)=deg(v)=h 



Ui,Vi 



u,v€V(G) i=l 
deg(u)=deg(v)=k 



E 

i=i 



E 



«,«ev(G) 

deg(u)=deg(v)=k 



Ei^o)i (fc) -i^,i)i (fc) 



i=i 



This completes the proof. ■ 

An advantage of Theorem 15.11 comparing to computing TWt{G) by the definition is that 
we do not need to compute distances, but only to count vertices in the classes. This theorem 
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can be considered as another instance of Klavzar cut method. For its general description and 
an overview of its applications in chemical graph theory see survey [19] . 

As an example, we obtain a closed expression for TW3 of the coronene/circumcoronene 
homologous series H k - In Figured! 2k — 1 horizontal elementary cuts of H k are presented. 
There exist two additional groups of 2k — 1 equivalent cuts, obtained by rotating the former 
group by ^ and — ^. The number of vertices of H k equals n k = 6k 2 , while there are exactly 
6k vertices of degree two. 




Figure 4: The Coronene / Circumcoronene H3. 

Using symmetry, the contribution of the elementary cut Ci is equal to the contribution 
of C2fc_j, i = 1, 2, . . . , k — 1. By induction it follows that for i = 1,2, ... ,k the number of 
vertices above cut Ci equals i{2k + i), while the number of vertices of degree 2 equals k + 2i. 
Therefore, 

k - 2i){<6k 2 - 6k - 2ki - i 2 + k + 2i) 

19k 
' Iff' 

Finally, we derive the fifth-order polynomial formula for the generalized terminal Wiener 
index of H k 

TW 3 (H k ) = -(k - l)k(2k - l){82k 2 - 82k - 19). 
5 

Using similar methods as in [TJ \T7\ [20] one can obtain the closed formulas for other 
chemical graphs (trees, benzenoid chains, phenylenes, ...) and design a linear algorithm for 
TWk of benzenoid systems. 

6 Concluding remarks 

The Wiener polarity index and the terminal Wiener index are very new molecular-structure 
descriptors and only a limited number of mathematical and chemical properties were estab- 
lished so far. In this paper we generalized these indices and open new perspectives for the 
future research. 

Another generalization of these indices may be the following 

W£(G) = \{(u,v) I d{u,v) < k, u,v G V}\ = W 1 (G) + W 2 (G) + ... + W k {G) 



fc-i 



-TW 3 {H k ) = (3k 2 -3k) 2 + 2^(2ki + 



8=1 



164/c 5 
15 



82k A 58k 3 



5fc2 
3 
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and 

TW* k (G) = Yl d(u,v)>TW 1 {G) + TW 2 (G) + . . . + TW k (G) . 

u,v€V(G) 
deg(u)<k, deg(v)<k 

It would be nice to study mathematical and algorithmic properties of these indices and 
report their chemical relevance. These indices are obtained from the famous Wiener index, 
which has many applications in chemistry, graph theory and computer science. 

Acknowledgement. This work was supported by Research Grant 144007 of Serbian 
Ministry of Science and Technological Development. 
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